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Abstract 

In the application of potential models, the use of the Dirac equation in central potentials remains 
of phenomenological interest. The associated set of decoupled second-order ordinary differential 
equations is here studied by exploiting the phase-integral technique, following the work of Froman 
and Froman that provides a powerful tool in ordinary quantum mechanics. For various choices of 
the scalar and vector parts of the potential, the phase-integral formulae are derived and discussed, 
jointly with formulae for the evaluation of Stokes and anti-Stokes lines. A criterion for choosing 
the base function in the phase- integral method is also obtained, and tested numerically. The case 
of scalar confinement is then found to be more tractable. 

PACS numbers: 03.65.Pm, 03.65. Sq, 12.39.Pn 
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I. INTRODUCTION 



Several problems of interest in theoretical physics lead eventually to the differential equa- 
tion 



( — +R(z)]4>{z) = 0, (f.l) 



\dz 2 

where R is a single- valued analytic function of the complex variable z. The form of (1.1) 
suggests looking for solutions expressed through a prefactor A(z) and a phase w(r), i.e. 



i)±{z) = A(z)e ±iw{z) . (1.2) 



The Wronskian of ip+( z ) an d i>-{ z ) is equal to —2iA 2 ^-, and on the other hand the Wron- 
skian of two linearly independent solutions of Eq. (1.1) is a constant. Thus, the prefactor 
A(z) reads as const, x— /===> an d one has [l| 

m = -jL=e^\ (1.3) 

where 

«;(*) = /%(CR, (1-4) 



the function w being the phase integral, while q is called the phase integrand. Moreover, 
upon insertion of the exact solution (1.3), (1.4) into Eq. (1.1), one finds that the phase 
integrand q(z) should satisfy the g-equation 

/ N „/ NN 3 d 2 1 R(Z) 

f(z, q(z), R(z)) = <f-— 2 q^ + -AJ - 1 = 0. (1.5) 

In practice, however, the task of finding exact solutions of Eq. (1.5) is rather difficult. The 
best one can do is often to determine a function Q that is an approximate solution of the 
g-equation (1.5), so that 

£ = f(z,Q(z),R(z))«l. (1.6) 

The approximate phase-integral method consists in finding approximate solutions of Eq. 
(1.1) with unspecified base function Q. A criterion for finding Q is that the function Eq 
defined in (1.6) should be much smaller than unity in the region of the complex- z plane 
relevant for the problem. However, this criterion does not determine the base function 
Q uniquely, the physicist has a whole set of basis functions Q at his disposal, and this 
arbitrariness can be exploited. 



On the other hand, along the years, many efforts have been devoted in the literature 
to the theoretical investigation of light fermions confined by a potential field [2|. In the 
phenomenological applications, when dealing with mesons consisting of a heavy quark and 
a light quark, one can imagine that the heavy quark is indeed very heavy and acts as a 
"classical" source that can be represented as a superposition of Coulomb-like plus linear 
potential, better known as Cornell potential j^]. The mass occurring in the Dirac equation 
is therefore the mass of the light quark. It is by now well known that, on using the Dirac 
equation, only Lorentz scalar confinement leads to normalizable stationary states, while in a 
suitable variant of the Dirac equation, called "no pair" , only Lorentz vector confinement has 
normal Regge behaviour. Hereafter we focus on the stationary Dirac equation for a quark of 
mass m in a Lorentz scalar potential Vs(r) and in the time component of a Lorentz vector 
potential VV(r), i.e. [4] 

dF _ n p | mc 2 + E + V s -V V c 
dr r he 

dG = !iG+m >-B + Vs + Vv 
dr r he 

where k = — I — 1 if j = I + \, K — I if j — I — \. In the resulting second-order equations, 
first derivatives can be removed by putting 




V# + me 2 + V s -Vv 

VE - me 2 -V s -Vy 




;i.9) 



Section II studies the second-order equations resulting from the radial Dirac equations 
(1.7) and (1.8), preparing the ground for the application of the phase-integral method. 
Section III describes various possible choices of basis function in the phase-integral method. 
Sec. IV arrives at a general criterion for choosing a suitable basis function Q. Sec. V 
performs a numerical analysis of the applicability of such a criterion. Sec. VI studies Stokes 
and anti-Stokes lines for the squared Dirac equation in a central potential, inspired by the 
choice of Q made in the simpler analysis of central potentials in ordinary quantum mechanics. 
Concluding remarks and open problems are presented in Sec. VII. 
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II. SECOND-ORDER EQUATIONS FROM THE RADIAL DIRAC EQUATION 



With the notation in the Introduction, our starting point is the following set of decoupled 
second-order equations obtained from the radial Dirac equation: 

d 2 



+ Rf(r) 




d 2 

dr' 2 





+ Ro(r 



where the "potential" terms read as J4] 




0. 



RAr) 



(E — V\ 



v 



{mc 2 + V s ) 2 k(k + 1; 



+ 



(hcf 



mc 2 + E + V s - V v 



K 



Vy-V^ 



r 2[mc 2 + E + V s - V v ] 
1 2 



mc 2 + E + V s - V, 



v 



(2.1) 



(2.2) 



R 9 (r) 



(E - V v ) 2 - [mc 2 + Vs) 2 «(k - 1) 



+ 



(he) 2 

Vy + V S 



E + Vs + Vy 



mc 2 - 



K 

r 



+ V 



2[mc 2 -E + Vs + Vy] 



Vy + Vl, 



,nr--E + Vs + Vy_ ' (2 " J>) 

Note that, for energies E < —mc 2 , the following difficulty arises: the relation between f(r) 
and F(r) in (1.9) becomes singular at the point r = r/ such that Vy(rf) —Vs(rf) = E + mc 2 . 
Thus, the effective potential R/{r) in (2.2) becomes infinite at r — > r/. The solutions become 
meaningless near the point r = rf because the phase integrals diverge. Similar remarks 
hold for g(r) and the effective potential in (2.3). However, this difficulty is purely 
formal because the original Dirac system (1.7) and (1.8) is not singular at the point r — Tf. 
A powerful JWKB analysis of the first-order Dirac system (1.7) and (1.8) can be found in 



[71 



Equations (2.1)-(2.3) suggest exploiting the known properties of the differential equation 
(1.1), which, as we said, is much studied in classical mathematical physics and ordinary 
quantum mechanics. The change of dependent and independent variable that preserves the 
form of (1.1) without first derivative is given by 

1 
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(2.4) 



w{z) = J Q(()dC (2.5) 

where the function Q is not specified for the time being but will be suitably chosen later. 
Upon defining 

-#- 1 + «- 3/2 |>- I/2 >> < 2 - 6 ' 

equation (1.1) can be expressed in the equivalent form 

d 2 



+ 1 + e 



(f(w) = 0. 



(2.7) 



Equation (2.7) is more convenient because it can be turned into a system of two linear 
differential equations of the first order. For this purpose, one assumes that the complex 
w-plane is cut in such a way that the functions appearing are all single-valued and (p can 
read as 

(p(w) = a x {wy w + a 2 {w)e- iw . (2.8) 

If we further impose that 



a[(w)e lw + a' 2 (w)e- lw = 0, 



(2.9) 



the first derivative of (p reduces to 



dip 
dw 



ia x e lw - ia 2 e~ lw , 



(2.10) 



and one obtains the desired system of two first-order ordinary differential equations, i.e. 

da± i 



dw 2 
da 2 i 



e <2i + a 2 e 



-liw 



dw 2 

Such a system can be written in matrix form as 

do 



e[a 2 + die 



2iu< 



dw 



M(w)a, 



having set 



2 ' 1 

M(w) = -e 

l \ _ 2iw _y 



2iw 



(2.11) 
(2.12) 

(2.13) 

(2.14) 



a\w) 



ai[w) 
a 2 (w) 



(2.15) 



At this stage, one can replace the differential equation (2.13) by the integral equation 

PW 

a(w) = a(wo) + / M (wi)a(wx)dwi, (2-16) 



J wo 

which can be solved by iteration, starting from the solution formula 

a(w) = F(w,wo)a(wo), (2-17) 

where 

pw rw t*wi 

F(w,w ) = 1+ / dw 1 M(w 1 )+ / dw 1 M(w 1 ) / dw 2 M(w 2 ) 

J WO J WO J WQ 

pw rw\ pW2 

+ / dw x M(w x ) l dw 2 M(w 2 ) / dw 3 M(w 3 ) + ... (2.18) 

J WQ J wo J WO 

Under the assumption that 

Y^\Mi3(w)\<m(w), (2.19) 
j 

where m(w) is a non-negative quantity, one finds that, in any region of the complex-w plane 
where the integral f m(wi)dwi is bounded, the series in (2.18) is absolutely and uniformly 
convergent. From (2.4), the original equation (1.1) is then solved by 

iw (2) — iuu (•s) 

^ = fll 4^fT + fl2 ^fT = a '(^i^ + a 2(^)/ 2 (2.20) 

vQ( z ) vQ( z ) 

where 

h{z) = 1. e lw{z \ Uz) = i e~ iw(z l (2.21) 

Our main source on this topic, ref. [8j, contains all details about useful approximate formulae 



for the F-matrix and many peculiar properties o 
should not be confused with the JWKB method 



the phase-integral approximation, which 



III. CHOICE OF THE BASE FUNCTION 



The function Q in Sec. II need not coincide, when squared up, with the function R in 
Eq. (1.1). A guiding principle in the choice of base function is as follows: first find the pole 
of higher order (if any) in R(z), and then choose Q(z) in such a way that it cancels exactly 
such a pole (see below). 
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A. Scalar confinement 

For example, the scalar confinement is achieved with the potentials 

V s = ar, V v = 0, (3.1) 

for which the "potential terms" Rf and R g in (2.2) and (2.3) reduce to 

^ E 2 — (mc 2 + ar) 2 k(k +1) a At 3 a 2 

^ (he) 2 r 2 (mc 2 + E + ar) r 4 (mc 2 + + ar) 2 ' 



2 



^ E 2 — (mc 2 + ar) 2 k(k — 1) ^ a At 3 a [3 3) 

9 (ac) 2 r 2 (mc 2 — E + ar) r 4 (mc 2 — E + ar) 2 



The experience gained in ordinary quantum mechanics suggests therefore choosing 

Q^ Rf+ (3.4) 

At (At — 1; 



Ql = R 9 + '^—T^- (3-5) 



B. Logarithmic potential 

More generally, however, bearing in mind that singularities in (2.2) and (2.3) might receive 
a further contribution from Vy or Vs if they were of logarithmic type, one can take 

V s = -\og(-), W = 0, (3.6) 



a \r 0/ 

bearing also in mind that only a scalar potential is able to confine a quark in the Dirac 
equation, and that a relativistic Qq system is indeed well described by the choice (3.6), as 
shown in 9]. The potential terms Rf and R g in (2.2) and (2.3) are then found to develop 
also a logarithmic singularity at r = 0, because the l'Hospital rule for taking limits implies 
that 

lim — — — = lim 5 = 00. 

r->o r 2 log(r) r->o r 2 log (r) 

We are then led to get rid of both the pole-like and logarithmic singularities of Q at r = 0, 
by defining 



«( «+i) + («+j) + j_ 

r 2 r 2 log(r/r ) 4r 2 log 2 (r/r ) 

At(At-l) 3 1 



Q 2 9 = R 9 + ' r 2 - r 2 V log ( r / ro ) + 4r^log 2 (r/r )' ( ' J) ' <S) 



Interestingly, we are suggesting a novel perspective on the logarithmic potential, arriving at 
it from the point of view of the singularity structure of the base function in the phase-integral 
method. 



C. A linear plus Coulomb-type potential 

One can also consider the Cornell potential [3( which is linear in the scalar part and of 
Coulomb-type in the vector part, i.e. 

V s = ar, W = (3.9) 
r 

As r — > 0, the centrifugal term in Rf (respectively R g ) is then found to receive further 

contributions with a second-order pole at the origin, so that we can remove such a singularity 
in Q by defining 

n 2_ R [* 2 - 1 - {b/ncf] 

Qf = H f + - 2 > 

Qg = R 9 + ~ 2 • 

However, the resulting integral (2.5) for the independent variable w is too complicated for 
analytic or numerical purposes. 



D. Analogy with central potentials in ordinary quantum mechanics 

It is therefore more convenient, in our relativistic problem, to fully exploit the arbi- 
trariness of the base function Q by defining it in such a way that it coincides with the 
form taken by Q in non-relativistic problems in a central potential. For example, for the 
Schrodinger equation in a central potential it is helpful to deal with a Q function of the 

n 

form [8J Q 2 = 1 + — . In our problem, both Rf in (2.2) and R g in (2.3) contain exactly, i.e. 
without making any expansion, the term — which is indeed of the form ^2 with 

(3,0, 



{he) 2 ' 

We thus look for 



Q 2 (r)=R f (r) + u f (r) = l + ^. (3.11) 



r 



S 



In this equation, the desired additional term can be obtained in exact form as 

In 

u,(r) = l + -/- J R/(r), 



where Rf leads to exact cancellation of the terms proportional to We then find, from 
(2.5) and (3.11) (see 




r j r 
2 V + \I 1 + 2 V 



(3.12) 



and, from (2.6), 



Rf -2 d 2 _i Rf _i d 2 i 

£f = Q 2 ~ 1 + Qf2 d^ Qf2 = Q}~ 1 ~ Qf2 d^ Q2f " 



(3.13) 



which yield, by virtue of (2.21), 



0/ (r) 



(3.14) 



where the functions aij and a2,/ can be obtained from (2.13)-(2.18), with e = Ef in (2.14). 
By following an analogous procedure, we find 



Q 2 g (r) = R g ( r )+u g (r) = l + ^ = Q 2 f (r, 



w g( r ) = w f( r ), 
R _a d 2 -i R —id 2 1 

£9 = of " 1 + Qs 2 2 = of " 1 " Qa 2 ' 



#0) = a i, fi 



3 iuj g (r) 



-iu>g(r) 



VQ 6 



+ a 2 ,, 



(3.15) 

(3.16) 
(3.17) 

(3.18) 



bearing in mind that R g ^ Rf = 
evaluation of a^ g and a 2jS . 

We should now recall that, by virtue of the identity 



e g 7^ Ef, and setting now £ = E g in (2.14) for the 



M{w 1 )M{w 2 )...M{w r) 



- j e(w 1 )e(w 2 )...e(w n ) [l - e^-"*)] [i - e -w(^-*a)] 



h _ g— 2i('Ui n _i-i«„)j 



-2iw n 



_ e 2iw 1 _ e 2i(w 1 -w„) 



(3.19) 
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i.e. |s 



the F-matrix in (2.17)-(2.18) can be expressed through a fairly simple series, i.e. [8] 

f w i 
Fn(w,w ) = 1 + / dw! -e{wx) 

J W * 



rW ^ rwi £ 

+ / dwi -e(iwi) / <iw 2 -£(^2) [1 - e 

Jwo * "'ion * 



'too ^ 

+ (3.20) 



i 



F 12 (u>,w ) = / dwi -e^Oe" 2 ™ 1 

'wo Z 



pw ■ pwi ■ 

\ d Wl -e( Wl ) / dw 2 -e(w 2 ) [1 - e -*(«*-«*>] 

Jwo * «/«)n ^ 



2i(uii-u) 2 )] ^—2iw2 



>w " Jwo 

+ (3.21) 

pw ■ 

F 21 (w,w ) = - / dw x -e{ Wl )e 2iwi 

Jwq ^ 

d Wl -eMe 21 ^ / dw 2 -e(w 2 ) [l - e" 2 ^-^)] 

Wo J wo 

+ (3.22) 

r % 

F 22 (w,w ) = 1 - / dwi -e(wi) 

Jwo ^ 

w ■ rw! ■ 

dwt -e( Wl ) / dw 2 -e(w 2 ) [1 - e~ 2l ^- W2) ] e 2i( - w ^ 

Wo J Wo 

+ ... . (3.23) 



IV. A GENERAL CRITERION FOR CHOOSING THE BASE FUNCTION 

We have also tried to find a base function Q by assuming its behaviour for small and 
large values of r, i.e. 

Q(r) = — + a 2 + a 3 r. 4.1 
r 

This base function can be analytically integrated, thus, in principle, we can obtain the phase 
integral according to (2.5). To fix the free parameter entering the previous expression we 
assume that the e parameter in (2.6) should vanish at small and large distances. However, 
this criterion does not ensure that e remains small throughout the whole range of values of 
r, and we have instead found regions where the resulting e is, regrettably, larger than 1, 
thus making our choices unsuitable. A general method is instead as follows. Since we have 



10 



to fulfill the condition (1.6) with e defined as in (2.6) and R = Rf or R g , we re-express 
in the form 

d 2 

' ' «Q 2 , 



R-Q 2 + Q l/2 ^-JQ- l/2 ) 

dz 2 



and define 



or, the other way around, 



bearing in mind that 



A = R-Q 2 , 
B = Q 1/2 ^(<T 1/2 ), 



^Q 1/2 |^(Q- 1/2 ), 
B = R-Q 2 , 



\A\ - \B\ < \A + B\ < \A\ + \B\. 



Moreover, we can always make the conventional choice according to which |^4| > \B\. 

When (4.3) and (4.4) hold, if both A and B are positive, the conditions (1.6) and 
yield 



i.e. 



R-Q 2 + Q 1/2 ^(Q- 1/2 )<Q 2 , 

dz 1 



R-2Q 2 + Q 1 / 2 ^(Q- 1 / 2 )<0. 
dz 1 



When (4.3) and (4.4) hold, if A > and B < 0, conditions (1.6) and (2.6) yield 



\R-Q 2 \ - 



<Q 2 , 



1.6) 
4.2) 

4.3) 
4.4) 

4.5) 
4.6) 

4.7) 

2.6) 
4.8) 

4.9) 



(4.10) 



which coincides with (4.9) because A = R — Q 2 > while B = —\B\ < 0. 

Nothing changes if instead (4.5) and (4.6) hold. For example, if A defined in (4.5) is 
positive and B defined in (4.6) is negative, one finds from (1.6) and (2.6) 



Q 1/2 ^- 2 (Q- 1/2 )-\R-Q 2 \<Q 2 , 

dz 1 



(4.11) 



which coincides with (4.9). Thus, in all possible cases, the family of as yet unknown base 
functions Q has to be chosen in such a way that the majorization (4.9) is always satisfied. 
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012345 012345 



r(GeV _1 ) r(GeV _1 ) 

FIG. 1: The left-hand side of eq. (14. 9b is plotted versus r for Sec. IIIA. The continuous (dashed) 
line corresponds to k = — 1 (k = 0). We have used m = 0.300 GeV, a = 0.308 GeV 2 and 
E = 1.9 GeV. 



V. NUMERICAL RESULTS ON Q 2 „ 

In this section we collect all numerical results regarding the choice of the squared base 
function Q 2 by following the considerations in the previous sections. First of all we work 
in the natural unit system (h = c = 1), and we plot in the figures 1-4 the left-hand side of 
eq. (14. 9p . The chosen range for r is the typical one for the heavy mesons phenomenology. 
The numerical values for the parameter are taken from the phenomenological analysis of 
the meson spectrum by using the Dirac equation [2]. In particular, we restrict ourselves to 
consider the numerical parameter for the charmed particles. Moreover, it should be observed 



that in [2] only the Cornell potential has been considered (cf. subsection IIII Cp . However, we 
use the same numerical values for parameters also in the case IIII A[ IIII CI and IIII Dl because 
the qualitative behaviour of the results does not depend strongly on the numerical values of 
the parameters. 

In figure [1] we have plotted the left-hand side of eq. (14. 9 j) for the (R, Q 2 ) = (Rf, Q 2 ) 
(left panel) and (R,Q 2 ) = (R g) Q 2 ) (right panel). The light quark mass, m = 0.300 GeV, 
a = 0.308 GeV 2 and E = 1.9 GeV in R f and R g (cfr eqs. flO)- (ET3]) ). The plots in figures 
20 are obtained by using the values collected in their captions. It should be observed that 
in figure [1] we have used a confining linear potential and for Q 2 the choice in section IIII Al 
The inequality in (14. 9 j) is satisfied for almost the whole physical range of r. 

In figure [2J the logarithmic potential has been considered (cf. section IIII Bp with ro = 



1 GeV -1 . Also in this case we do not have direct phenomenological information on the values 
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012345 012345 
r(GeV _1 ) r(GeV _1 ) 

FIG. 2: The left-hand side of eq. (|4.9f) is plotted versus r for Sec. IIIB. The continuous (dashed) 



line corresponds to k = —1 (k = 0). We have used m = 0.300 GeV, a = 1/V0.308 GeV" 1 , 
E = 1.9 GeV and ro = 1 GeV. The left (right) panel corresponds to the case of / (g) in eq. (|2.ip . 
Note the range of r. Moreover, it should be noticed that, for r G [0, 1], the inequality in eq. (14. 9h 
is strongly violated. 




FIG. 3: The same as figure 1. Moreover, b = —0.579. 



of the parameters. Smaller values for ro are responsible for the violation of the inequality 

(USD. 

In figures [3] and H] the Cornell potential is considered. In these figures the values of the 
parameters are taken, as already said, from the phenomenological analysis. In fig. [3] the 
inequality is violated for Q 2 g in the whole range of r. While the case inspired by ordinary 
quantum mechanics (cf fig. H]) violates the inequality in the region of small r. 
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r (GeV -1 ) 



FIG. 4: The same as figure 3. Here Q 2 g are chosen as in section fill Dl 



VI. STOKES AND ANTI-STOKES LINES 



In the application of the phase-integral method to Eq. (1.1), a concept of particular 
relevance is the one of Stokes and anti-Stokes lines. By definition, the differential dw = 
q(z)dz (see (1.4)) is purely imaginary along a Stokes line, and real along an anti-Stokes line. 
Thus, the Stokes lines are lines along which the absolute value of e tw ( z ' increases or decreases 
most rapidly, while the anti-Stokes lines are level lines for constant absolute values of e lw ^ 

Q. 

For example, for the case studied in our subsection 3.D one can evaluate at complex 
r = Ae %e the phase integral (3.12). One then finds, after repeated application of the Gauss 
representation of complex numbers, and upon defining 



A = l + Acos(0), B = ^sin(9), 



6 = arctan 



Asm{6) 



(277 + Acos(0)) 
B 2 



A 9 , 



cos • 



sin 



2' 


2' 



sm 



ip = arctan -p^- 



#-cos|+(A 2 + 5 2 ) 4 cos 



(6.1) 
(6.2) 
(6.3) 
(6.4) 

(6.5) 
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the following split of Wf(r) into real and imaginary part: 

'e-e 



Re Wf = 2rj 




( A 2 + B 2 Y cos { U -^-\ + ±k>g(n 2 + 



Imwf = 2r) 




(6.6) 



(6.7) 



From what we said before, along an anti-Stokes line, dwf is real, and hence Imwf is constant. 
We thus find from Eqs. (6.5) and (6.7) the transcendental equation 



/ #sin|+U 2 + £ 



arctan-^- 



sm 2 



^-cos^+fA 2 + 5 2 Vcos 




A (a 2 + b 2 " 



' 9 + 9 . 
sin | — - — I = const. (6.8) 



2t) 

Moreover, since dwf is purely imaginary along a Stokes line, we are led to consider the 
equation 

Re Wf = constant. 

This becomes, from (6.6), the transcendental equation 




+ log(« 2 + (3 7 



const. 



(6.9) 



In general, we cannot give analytical solutions to the equations (16.81) and (16.91) . However, 
the fact that, for reasonable values of the parameters, solutions to such equations exist is 
crucial. In this respect, in figs. [5] and [6] we show that, for r\ = (0.5, 1, 2) and A = 3, they can 
be solved for a constant value and for zero, respectively. In particular, eq. ( 16. 8ft has either 
zero or six roots depending on the choice of the value for the constant, unlike the case of eq. 
( 16. 9p . where at most three zeros can be found depending on the constant. 

Following what we say at the beginning of this section, the absolute value of e iWf increases 
of decreases along the Stokes lines while it remains constant along anti-Stokes lines. Figure 
M displays this behaviour in a neat way. 



VII. CONCLUDING REMARKS AND OPEN PROBLEMS 



Second-order equations for relativistic systems have been investigated a 
including the work in 10], and supersymmetric extensions considered in 



ong many years, 
111 ]. In ordinary 



15 



3r 

2 - 




FIG. 5: Here we plot the left-hand side of eq. (|6.8p minus 2, which is the value of the constant in 
the same equation, versus 6. The curves have been obtained for 7] = (0.5, 1, 2) (dashed, continuous, 
long dashed) lines (cf eq. (|3.10p ) and A = 3. This figure shows that solutions to eq. (|6.8p exist. 








FIG. 6: Here we plot the left-hand side of eq. (|6.9p versus 9. The plot is obtained for 77 = (0.5, 1, 2) 
(dashed, continuous, long dashed) lines (cf eq. (|3.10p ) and A = 3. As in the case of fig. solutions 
to eq. (|6.9p exist. 
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aantum mechanics, the most powerful choice of base function is the Langer choice 
, Q, but the peculiar technical difficulties of the effective potentials (2.2) and (2.3) for 
the Dirac equation cannot be solved in the same way, and one has rather to resort to the 
JWKB method along the lines in [zj]. It was here our intention to investigate potentialities 
and limits of the phase-integral method, which actually differs from JWKB methods [l|. 
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FIG. 7: We plot the left-hand side of eq. (|&8j) versus 6 and A for 77 = 2 GeV. 




< 0.2 - 



FIG. 8: Here we plot \e lw f\ versus 9. The curves have been obtained for rj = (1,2) (continuous, 
long dashed) lines (cf eq. (|3.10|i ) and A = 3. 



Our results are of qualitative nature, while we fail to obtain bound-state energies from the 
integrals in sections 2 and 3. At a deeper level, the problem arises of solving coupled systems 
of first-order ordinary differential equations which, when decoupled, give rise to a pair of 
equations of the form (1.1). The phase-integral method, originally developed for second- 
order equations of the form (1.1), should have implications for the solutions of the original 
first-order system as well. This expectation should be made precise, and its relation with 
the JWKB method should be elucidated. 

Although the decoupled second-order equations obtained from the radial Dirac equation 
are formally analogous to the second-order equations to which the phase-integral method 
can be applied, the actual implementation is much harder because the "potential" terms Rf 
and R g therein contain complicated denominators built from the potentials V$ and Vy in 
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the radial Dirac equation [4|, [15| . This implies that the actual choice of base function Q is 
a difficult problem. In Sec. Ill we have described some possible choices of Q, and in Sec. 
IV we have arrived at the majorization (4.9) to select Q, tested numerically in Sec. V. The 
analysis of (4.9) for the Cornell potential shows that an appropriate basis function can be 
found for the case k = —1 (see fig. 3). Moreover, for the logarithmic potential the plots 
displayed in fig. 2 show that (4.9) is not fulfilled in the whole range of r. The investigation of 
Stokes and anti-Stokes lines in Sec. VI is also, as far as we know, original in our context. It 
remains to be seen, however, whether such lines can be of direct phenomenological interest. 

The work in Ref. [4|], despite being devoted to the amplitude-phase method, did not 
investigate our same technical issues. Thus, no obvious comparison can be made. The years 
to come will hopefully tell us whether choices of Q satisfying (4.9) exist for which the F- 
matrix in (2.18) can be actually evaluated. In the affirmative case, one would gain conclusive 
evidence in favour of the superiority of the phase-integral method. In the negative case, one 
would instead gain a better understanding of the boundaries to our knowledge. 
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